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We propose an approach which allows to construct and use a potential function written in terms of 
an angle variable to describe interacting spin systems. We show how this can be implemented in the 
Lipkin-Meshkov-Glick, here considered a paradigmatic spin model. It is shown how some features 
of the energy gap can be interpreted in terms of a spin tunneling. A discrete Wigner function 
is constructed for a symmetric combination of two states of the model and its time evolution is 
obtained. The physical information extracted from that function reinforces our description of phase 
oscillations in a potential. 



Quantum tunneling is a very peculiar quantum phenomenon, and has long been studied in a variety of contexts. 
Due to different aspects of the process, coherent or incoherent tunneling have also been distinguished, according to 
a possible periodicity of the motion of the particle trapped in a potential. In the first case, the motion leads to a 
characteristic oscillation frequency related to the tunneling splitting of the ground state energy, whereas in the second 
the particle will be find outside the confining potential region with a mean scape rate. In this last case, semiclassical 
methods are usually employed in order to calculate the mean scape rate and tunneling probabilities, although other 
methods are sometimes of great interest [1]. 

In what concerns the spin degree of freedom, it has been shown that quantum tunneling also appears in several 
physical situations such as magnetic moments in spin glass systems, single magnetic impurity in a crystal field, a single 
domain, a ferromagnetic, or antiferromagnetic, domain wall, and so on [2]. In these cases, transitions between states 
associated to energetically equivalent configurations exist, even at low enough temperatures to exclude thermally 
assisted transitions. Those transitions were then suggested to be due to quantum tunneling. 

From the theoretical point of view, the general approach to the spin tunneling process was first treated by van 
Hcmmen and Siito [3] , and by Enz and Schilling [4] . Further developments came along the approaches which make 
use of a generalization of the WKB method adapted to spin systems [3,5-7], and, alternatively, by using the Feynman's 
path integral formulation of quantum mechanics [2,8]. In general, the starting point of any of the two approaches 
is a Hamiltonian that describes a spin degree of freedom which interacts with an external/internal field, for instance, 



associated to the manganese acetate complex [2], being these Hamiltonians written in terms of the standard spin 
operators obeying the usual commutation relations 



In the present paper we propose an alternative useful approach to the description of the problem of quantum spin 
tunneling. The aim of this approach is to obtain a potential function, written in terms of an angle variable, which can 
account for the main properties of the spin system; in particular we intend to show that this potential can describe 
the proper behavior of the spin system when it undergoes a spin tunneling phenomenon. 



I. INTRODUCTION 
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In order to test our proposal we applied our approach to the Lipkin-Meshkov-Glick (LMG) model [9], which is 
defined through a many-spins Hamiltonian containing a spin-spin interaction term besides a one-spin term, i.e., a J z 
contribution. In fact, that Hamiltonian has been originally proposed to represent a system of iV interacting fcrmions 
in such a form that, due to the algebraic structure involved, it is a soluble model. It has been also discussed as a bona 
fide spin model mainly because of its phase transition in the TV — > oo limit [10,11]. Here we will be interested in its 
behavior for finite values N , mainly those values which can characterize mesoscopic systems. 

The starting point in the procedure to construct that potential is the Generator Coordinate Method (GCM) 
proposed long time ago in order to treat many-body systems [12] and also adapted to spin systems [13,14]. To carry 
out the calculations we choose the su(2) coherent states as our generator states and use the LMG Hamiltonian to 
describe the spin system. As we get the GCM energy and overlap kernels we construct a Hilbcrt space which stands 
for a new suitable kinematical ground to describe the spin system collective states. This can be achieved through a 
procedure put forth several years ago in order to avoid undesirable diverging solutions of the original GCM solutions 
[14,15]. Through a projection of the GCM energy kernel onto this new state space an equation for the spin system 
is obtained which presents an energy kernel that is a nonlocal function of the angular momentum. From this new 
energy kernel we can finally extract the potential function in terms of an angle variable [16]. 

The behavior of the LMG model gap curve as a function of the two-spin term interaction strength is then discussed 
in connection with the associated potential function, and we show that some aspects of the energy gap curve can be 
put in relation to potential curves which exhibit the characteristic behavior of spin tunneling. 

As a further test of the description of this spin system based on a potential function, we also constructed the 
Wigner function associated to a combination of the two lowest states of the model for some particular values of the 
strength parameter. For each case its time evolution was then obtained and the angle distribution probability function 
extracted. We then show that the results are also compatible with the potential function description. 

This paper is organized as follows. In section II we recall the LMG model as a spin system model and discuss some 
of its properties, mainly those relevant to the discussion of the spin tunneling. The extraction of the potential for the 
spin system is discussed and obtained in section III. In section IV we present the Wigner function for the two lowest 
energy states for different values of the interaction strength parameter and its time evolution is obtained. Finally 
section V is devoted to the final considerations. 



The Lipkin-Meshkov-Glick (LMG) soluble model [9] was originally proposed in order to test mean-field approxi- 
mations in many body systems, and was introduced as a Hamiltonian containing a one-body term plus a two-body 
interaction term in such a way that a set of N p (here considered an even number) interacting fermions is distributed 
in two levels separated by an energy e, being each one of them Ap-fold degenerate. The degenerated states within 
each level are labeled by a quantum number r, which assumes values between 1 and N p , and also by another quantum 
number a which is equal to +1 or —1 designating the upper or the lower level. The Hamiltonian of the model then 
reads 



Although the model was originally tailored for the treatment of fermion particles, it was soon verified that it could 
be put in a more simple form if we realize that a subjacent algebraic structure can be identified within the operators 
constituting the Hamiltonian. If we introduce the so-called quasi-spin operators defined by 



II. THE LIPKIN-MESHKOV-GLICK SOLUBLE MODEL 



A. The Lipkin-Meshkov-Glick Hamiltonian 
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we see that the LMG Hamiltonian can then be written as 

H = eJ z + ^(Jl + Jl), 

which exhibits a typical spin character. It is also direct to verify that those new quasi-spin operators indeed obey the 
standard su(2) commutation relations, viz., 

[J+,J_] = 2J Z , 
[J Z ,J±] = ±J±. 

In this way, if we forget for a while its original proposal, the LMG Hamiltonian can also be studied as a legitimate 
spin model of its own. The correspondence between the two operator languages is then direct: In terms of the fermions 
it is clear that the J z term gives half the difference of the number of particles in the upper and lower levels, while the 
second term, involving Jr. and J_ , is associated to the interaction between a pair of particles in the same energy level, 
and scatters this pair to the other level while preserving the quantum number q of each particle. In what follows we 
will treat the LMG Hamiltonian as a bona fide spin model. 

It is useful to rewrite that Hamiltonian in such a way to scale the interaction term to the particle number while 
measuring the energy in terms of e, i.e., 

with x = N p V/e. 

In what refers to the solutions of this model, we see that from general symmetry considerations some results can 
be advanced. First we observe that, with the Hamiltonian written in terms of these quasi-spin operators, and since 
[Hl,J 2 ] = 0, it can be diagonalized within each multiplet labeled by the eigenvalues of J 2 and J z , what accounts 
for the soluble character of the associated quantum problem. Each multiplet is 2 J + 1-dimensional, and we also note 

that the ground state belongs to the multiplet J = N p /2 = max J Zl being the J 2 and J z eigenvalues ^ + lj and 

— -^i respectively. Furthermore, we immediately see that the model also presents two discrete conserved quantities 
which reflect in the properties of the solutions. The simplest operator commuting with the Hamiltonian, therefore 
giving a constant of motion, is the parity operator 

n = cxp (i-K J z ) , 

indicating that the Hamiltonian matrix, in the J z representation, breaks into two disjoint blocks involving only even 
and odd eigenvalues of J z respectively. A second interesting property of the energy spectrum is revealed by the fact 
that the Hamiltonian anticommutes with the operator 

N p 

R = cxp (j-T^Jz^j CX P (™Jy) = I — m ) (— l) Wp+m (to I . 

m— — N p 

This operator is associated to a rotation of the angular momentum quantization frame by the Euler angles (— ^ , tt, 0) , 
thus transforming Hl into —Hl- As a result of this anticommutation property we observe that if | Ej) is an energy 
eigenstate with eigenvalue Ej, then the state R \ Ej) is also an eigenstate of the Hamiltonian with eigenvalue — Ej. 
This symmetry property of the Hamiltonian gives rise to an energy spectrum that is symmetric about zero. 

B. Some features of the energy gap 

In what concerns the exact solutions of the LMG Hamiltonian, they have long been discussed and found (some of 
them algebraically) for different -/V p 's by Lipkin [9]. They can also be numerically obtained by just diagonalizing the 
matrix associated to the Hamiltonian in the J z basis. Once the energy spectrum is obtained it is a simple task to find 
the energy gap, A = (Ei — E a ) /e, which is a very important quantity characterizing the behavior of this model as 
a function of the interaction parameter \- One of the main features of this model is that it presents a second order 
phase transition when N p — > oo, as already widely discussed in the literature, as can be seen from the energy gap 
curve in the N p — > oo limit which presents a typical discontinuity at x = 1, namely 
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lim 



(i-x 2 ) 1/2 Ixl < 1 

N p ^oo [0 \ X \ > 1 

Although this result has a particular interest in the study of systems presenting phase transitions [11], here instead, 
we will be concerned with the study of the model for finite, medium range values of N p since it can then stands for a 
mock mesoscopic spin system. 

In this context the main features of the LMG model can be directly extracted from the energy spectrum. As an 
illustration of this, Figure 1 shows a typical energy spectrum for N p — 10 as a function of the parameter \, the 
positive part of which is obtained just by reflecting the spectrum about zero, as mentioned before. As a direct result 
the energy gap curve can be obtained and is shown in Figure 2. It is clearly seem that the energy gap decreases as the 
parameter \ increases in such a way that we can infer the two-body correlation effects introduced by the interaction 
term of the Hamiltonian; in other words, as the two-body correlation strength increases the energy gap decreases. 
There are, however, as suggested by the shape of that curve, particular values of the strength parameter \ which 
indicate a change in the effects of the two-body correlations, namely we can extract the values of \ for which there 
occur those changes by just evaluating the first few derivatives of the energy gap curve. In fact, by obtaining the first 
numerical derivatives of the curve shown in Figure 2 we observe that there are in particular two values of \i as can 
be seen in Figure 3, that indicate those changes. These two values of \ separate the energy gap curve in three regions 
- as also shown in Figure 2. The first region can be understood as that one in which the mean field - controlled by 
the J z term - dominates, in contrast to the third region where the two-body correlations are the main responsible for 
the behavior of the spectra. The second region is characterized by the competition of both contributions. 

In what follows we will present a description based on a collective potential extracted from a full quantum mechanical 
formalism that can be used to visualize the main aspects of the regions described above. 



III. THE LMG COLLECTIVE POTENTIAL 



As already discussed in Ref. [16] a procedure can be developed for the LMG model that permits to treat the quasi- 
spin system as if it behaves as a particle system. In other words, we can associate a potential function of an angle 
variable defined in the interval — ir < <j> < n to describe the spins system interaction. 

For the LMG model the whole procedure can be made explicit by firstly assuming the su(2) coherent states as the 
generator coordinate method [12,14] variational space, as has long been proposed [13], namely 

I «) = cos*' (|) cxp (tan|j+) | 0), 

where | 0) is the vacuum of the angular momentum operator. The generator coordinate method then requires the 
calculation of an overlap kernel which, in the present case, is of the form 




Furthermore, an energy kernel must also be calculated which, in terms of the new variables ip = iH±£i and 9 = a — a, 
reads 

(a \H L \a) =H{ip,6) = 



N, 



p cos N - I - 



COS (if) cos 1 ( - 



cos 



-) (1 + sin 2 fo,))-! 



(2) 



A new discrete orthonormalized finite-dimensional representation, instead of this nonorthogonal continuous one, 
can be achieved [15] if we find a transformation that diagonalizes the overlap kernel, Equation (1). This can be 
directly accomplished by a Fourier transformation 



/" 

J — - 



(a | a) u m (a) da' = \ m u m (a) , 

r 

where the eigenfunctions of the overlap kernel are 
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, . cxp (ima) 
u m (a) = 



V2tt 

with the corresponding eigenvalues 



N (a\ ,. , , 2irN p \ 
I cos p — cxp (ima) da = -, r - i - 7 r— , 



with — < m < Therefore, the diagonalization process generates a set of N = N. p + 1 orthonormalized states 
defining a basis in the ground state multiplet we discussed before. The next step consists in projecting the LMG 
energy kernel, Equation (2), onto this angular momentum basis states, which leads to 



/7T /*7T 
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H(m,m) = / / exp [iip (to'— m)] exp 



-i— (m + to) 



H (<p, 0) dipd6 

\f\nKi 27T 



The diagonalization of this matrix gives the exact spectrum of the LMG Hamiltonian, as expected. 

Since the angular momentum space is characterized by a finite set of states, we can obtain a new discrete repre- 
sentation of the LMG Hamiltonian, viz., the discrete angle representation associated to the angular momentum one 
through a discrete Fourier transform, namely, 

N-l N-l 

exp (imdk) » exp (—irriOk) 



where now Ok labels the discrete angle variable, and 6k — ^ff-- 

The hnal entangled matrix of the energy kernel in this discrete angle representation does not allows us a direct 
comparison with Equation (2), written in terms of the continuous angle. However, the finite-dimensional character of 
this description can be exploited to allow for a discrete phase space description of this kind of degree of freedom [19] 
in which an angle-angular momentum pair labels each point of the space. One virtue of this phase space picture is 
that it further allows for the extraction of a potential function, described in terms of the angle variable, much in the 
same spirit as for particle systems, which can then describe the behavior of the spin system. This extraction can be 
implemented if we realize that the zeroth moment of this matrix, with respect to the variable 9k — Ok; represents a 
potential function in the discrete angular variable (Ok + Ok) /2. To this end, let us define 

9k — Ok' = u, 

and 

k + Ok- = 2<j). 

Due to the periodicity of the functions involved, we must only consider the interval for the summation over k to run 
as — < k < Furthermore, in order to emphasize the angular character of this discrete variable we will 

consider that range to be defined as ~^^-tt < u < ^r-7r, being the steps of this variable Am = 
In this form, the zeroth moment will be calculated as 



N — l 

M oW = E £ E / / 

.„-_N-±„ - N p J-7TJ-7T 



exp [i ((f) — ip) (m — m)} 



x cxp 



i(u-0) 



27r\/A m A 

(m'+ to) 



H (ip, 0) dipdO. (3) 



Defining the new variables r = rri— to and s — ( TO '+ m ) ; an( j us ing the general result 

^ TP -1 Np+r Np Np-r N p 
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__ N P _>_ N P r=-N p s=-Np-r r=l s=r-N p r=0,s = -N p 
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where the summations over s are restricted to run only over even(odd) values depending if r is even(odd), we can 
perform all the integrals and summations. Therefore, Equation (3) gives us the discrete potential function associated 
to the Lipkin model in terms of the angular variable <j). 

It has been also shown [16] that the large N p limit can be obtained from this approach which coincides with the 
result by Holzwarth [13]. The form for the potential, when N p » 1, is given by 

V (0) = M (0) = -fcll cos (0) - x(A ^ +3) sin 2 (0) . 

It is interesting to realize that the Lipkin model attains a classical limit when N p — > oo, in the sense that quantum 
dynamics becomes the classical dynamics [11]. Although this result guides us in the study of the behavior of large 
systems, for small or intermediate values of N p , when the two-body correlations are as important as the mean field 
term, and must be taken into account, the potential function must be calculated directly from Equation (3) which 
keeps all information regarding the spin system. 

Now we can look at the potential functions associated to the LMG model calculated for those special values of the 
interaction parameter \, obtained by the first three derivatives of the gap curve, as mentioned before. In Figure 4 we 
show the various curves associated to the potential for different values of \ f° r N p — 10; it is also shown in each case 
the ground state and the first excited state energies, which give the energy gap. It is immediately seen that for \ = 0, 
when the spectrum comes from a pure J z term, the potential presents a single minimum at (f> = 0. The bottom of 
the potential then becomes more and more shallow as x increases leading to a narrowing of the energy gap, i.e., the 
energy levels approach each other due to the uncertainty principle since the potential is then more wide in the angle 
variable. This behavior characterizes the first region in the gap curve. Now, the potential changes from a minimum 
to a maximum at <j> = when \ — 0-85 for N p — 10. Soon after the maximum at the origin appears, more specifically 
around \ ~ 1.2, the first particular value of the interaction, a more pronounced barrier starts to show up at the origin, 
and the wave functions present reflection effects so that the energy gap changes its behavior; for this value of x starts 
the second region. We then observe that for the second value of x, namely \ ~ 1.8, the second energy level touches 
the top of the barrier, and tunneling begins to play an essential role in defining the energy gap. Therefore, the second 
region ends while the third region starts which consists essentially of spin tunneling through a barrier. In Figure 4 
we also show the plot for \ — 2.5, a paradigmatic case of the third region, where tunneling fully takes place. 

Thus, we see that a potential function written in terms of an angle variable can be constructed for the LMG 
quasi-spin model in such a way that its description of the system is consistent with the exact results, obtained from 
the direct diagonalization of the energy matrix in the angular momentum representation, and can then account for 
a natural description of the behavior of the energy levels as the interaction strength changes. Furthermore, this 
potential function description may then stands for a natural tool to treat spin tunneling processes. 



IV. THE TIME EVOLUTION OF THE WIGNER FUNCTION 



Another approach to the LMG model energy gap curve comes from the study of the time evolution of the Wigner 
function associated to the lowest energy levels. Since the finite-dimensional phase spaces basic ideas have already 
been presented in the past [17,19], we briefly recall how this scheme can be implemented. 

First of all we identify the finite set of states of the angular momentum which labels the LMG model ground state 



multiplet \^f-,m) as the eigenstates of the Schwinger [18] unitary operator U, i.e. 



U | -f,m) = cxp I — m I | 



(4) 



in such a way that we can assume | ^-,m) = \ u rn ). 

Since we know the LMG quasi-spin Hamiltonian, the mapped Liouville operator which governs the time evolution 
of the Wigner function on the discrete phase space - labelled by the angular momentum and angle pair (m, n) - can 
be directly calculated by means of the mapping 

C(u,v,r,s) = ±Tr{tf(u,v) [H, — ]}, 



where 



G (m, n) — > — ^— exp 

3,1 



N 



in(3 (j, l;N)-^- (mj + nl) + i^jl 
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is the operator basis proposed in [19], U and V are the Schwinger unitary cyclic shifting operators [18], and the phase 
(3 (j, I; N) performs all the mod N arithmetic involved in the given operator basis calculations. 
It has been shown [17] that the discrete mapped Liouvillian reads 

£ (u, v, r, s) = 2i ^ ^2 k sin ^( bc ~ ad ) cxp [i?r$ (a, 6, c, d; A^)] 

m,n a,b,c,d 

exp | [a (u — to) + b (v — n) + c (u — r) + d (r — s)] 
where h (to, n) is the discrete phase space mapped form of the LMG quasi-spin Hamiltonian, namely 

Mm, n) = m + X ^ + to) + m + l) - to) -m+l) cos | 2 n, 

and $ (a, 6, c, rf; iV) = -(3 (^a + c+ ^f-,b + d+ ^f-;N^J . 

We also recall that the von Neumann-Liouville time evolution equation for the density operator, for time independent 
Hamiltonians, 

is represented in the discrete phase space by the mapped expression 

d 

{ qIPw (u, v ; t) = ^2 C ( M ' v > r ' s ) P™ ( r ' s; *) ' ( 5 ) 

r,s 

where now p w (u, v;t) is the Wigner function associated to the state of the system we are interested in [19] . A solution 
to Equation (5) in the form of a series 

p(t) = p (t ) + (-*) (t - to) [H, p (to)} + 1 (-if (t - to) 2 [H, [H, p(to)}} + ... 

has its discrete phase space mapped expression written as [16] 

Pw (u, v;t) = Y J {CJC + H) (* - *o) £ («, v, r, s) + 



(-if (t - tof ^2 ( u ' v > x ' v) L ( x > f> r ' s ) + ••• } ( r > s ; f o) • 



(6) 



Therefore, since we are given the Liouvillian for the LMG model, we can obtain the Wigner function at any instant 
of time t by using the series representing the iterated application of the discrete mapped Liouville operator. It is 
also to be noted that this time evolution is then carried out as a linear process of composition of sums of products 
of the arrays characterizing the Liouvillian and the Wigner functions respectively, over all the sites constituting the 
associated discrete phase space of the model. From an operational point of view, we have to construct the arrays 
defining the Liouvillian for the LMG model and the Wigner function. Once the Wigner function is given at time to, 
the calculations with the series, Equation (6), will give us the propagated Wigner function at any instant of time. 

Furthermore, as has been already shown, the probability distribution for the angular momentum and for the angle 
variables, as a function of time, can be obtained directly from the Wigner function by means of a trace operation, 

N p 
2 

L(m;t) = ^2 Pw(m,n;t), 

N p 

n= f- 

and 

Np 
2 

®(n;t)= Pw(m,n;t) (7) 

Np 

m= £■ 

respectively. It must be noted that the angle variable is also discrete being characterized by 9 n — j^n. 
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A. Wigner functions associated to some of the LMG quasi-spin model states 



Since the general mapping expression 



(m, n) = -^Tr [G f (to, n) 0] 



holds for any operator O acting on the states of the corresponding state space, we will be concerned with the density 
operator p, in particular for pure states 

P = I i>}{ip \, 

when these states can be expressed in terms of the basis states of the unitary U operator, Equation (4), which are 
also eigenstates of the J z operator, namely, 



ip) = E Cfe | u k ). 



For these cases it is simple to show [19] that the Wigner function reads 



p w (to, n) 



- t 

N 2 



exp 



2ni . . 

— (mj + nl) 



I {3,1), 



(8) 



where 



fu,i)= E 



°p ex p 



2ni 



and the index [p + 1} just stands for the value of the sum p + I that are cyclically restricted to the range of the labels, 
e.g., p max + 1 = Pmin- Furthermore, since we have chosen the ranges of the labels j, I and p to be symmetric about the 
origin, we will obtain real- valued Wigner functions, and due to this choice it must be also noted that, in this context, 
the mod N phase /3(m, n; N) present in the operator basis will not contribute to the above expression. The advantage 
of having Wigner functions in that form is that if we diagonalize the LMG quasi-spin Hamiltonian in the J z basis we 
can obtain directly the coefficients Cfe in order to build a Wigner function using Equation (8). In particular, we can 
construct the symmetric and antisymmetric combinations 



and 



4 I u l) 



4) 



from the i-th and j-th eigenstates of the LMG quasi-spin model respectively. These states will be of great interest 
specially when dealing with the ground and the first excited state since then we know that the energy gap curve as a 
function of \ indicates the existence of three regions of interest for finite N p . From those two states we immediately 
obtain, by means of the direct use of Equation (8), the two corresponding Wigner functions, i.e., p s w (m,n) and 
p° (to, n) respectively. 

To carry out the time evolution of the Wigner function for any one of these possibilities, we have only to apply 
Equation (6). Furthermore, it is also direct to see that we can extract the superposition probability of the t = state, 
represented by the corresponding Wigner function (symmetric/antisymmetric), with a time evolved one. This will 
indicate us the existence of periods of oscillation between symmetric/antisymmetric states for specified values of the 
interaction parameter \- To obtain this probability as a function of time we have to realize that 
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p *f (*) = ^2Pw ( m > n ) pL {™, n; t) , (9) 

so that we have only to carry out the summations over the whole discrete phase space of the product of the Wigner 
functions. 

All the numerical calculations were carried out for 1200 steps of time, each one being At = 0.1, and the family of 
curves corresponding to the symmetric case is presented in Figure 5. In all of them we started with the symmetric 
Wigner function at t = 0, and calculated the probability of obtaining the symmetric Wigner function again at any 
time t. It is clear that for x — 0.0, when the Hamiltonian is a pure J z term, and the symmetric state only involves 
the ground state and its closest neighbor with no coupling, the corresponding Wigner function will oscillate back and 
forth with period 2tt, since the energy gap, A, is exactly 1 (it must be recalled that we are measuring energies in units 
of e). This corresponds simply to a free oscillation, since the spectrum is equally spaced, there are no coupling term 
in the Hamiltonian, and no tunneling at all. We then observe that, as soon as the interaction strength \ begins to 
increase, the oscillations change their frequencies in such a form that the basic relation 

.(x) - ^> 

is verified, as is expected for the particular case of a mixture of two energy eigenstates of the Hamiltonian. 

As a by-product of the time evolution of the Wigner function, we can also obtain the mean value of the phase 
variable, <j>, of the system state at any instant of time by just considering the average of the phase distribution 
function, Equation (7), with the corresponding Wigner function, in the same form as that prescribed in Equation (9). 
The results for the symmetric case are presented in Figure 6. From those plots it can be clearly seen that while the 
Wigner function oscillates between the symmetric and antisymmetric configurations periodically in time, the mean 
phase bounces back and forth in the potential. As \ increases and the barrier at the origin begins to play a role in the 
time evolution, the oscillation frequency in the potential diminishes and the Wigner function spends more and more 
time in one of the two minima, until it tunnels to the other one. For sufficiently high values of \ the frequency tends 
to zero (since A goes to zero), and the angle distribution keeps trapped into one of the two minima of the potential 
therefore stopping the tunneling process. 



V. CONCLUSIONS 



In this paper we have shown that a potential function description for the quasi-spin Lipkin-Meshkov-Glick model 
can be constructed that can account for all properties a typical spin system presents. Although originally based in the 
treatment of many particle systems, this description has been partially developed and exploited in the past, and shown 
here to be of use for spin systems as well. In fact, we have explored the LMG model for finite intermediate number of 
spins in order to simulate mesoscopic spin systems. As such, we have shown how to describe some interesting features 
of the energy gap curve exhibited in that model for N — 10 in terms of the behavior of the spin system phase state in 
a potential. In particular this description allows us to interpret the appearance of three distinct regions in the energy 
gap curve as the change of behavior of that phase state when the form of the potential changes as the interaction 
parameter \ increases. When \ exceeds a certain value, we are led to a typical tunneling process in the proposed 
potential that can then be fully described in terms of a phase variable. In this sense the proposed potential is an 
alternative useful tool to handle spin tunneling processes. 

We also have shown that a discrete Wigner function associated to a symmetric (or antisymmetric) combination 
of the ground and first excited states of the LMG model can be constructed and its time evolution obtained when 
considering the proposed potential. The time behavior of the Wigner function corroborates our description of the 
LMG quasi-spin model. The phase mean value can also be directly obtained from the Wigner function and clearly 
shows the oscillations of the phase state in the potential. When \ greatly increases, leading to the presence of two 
pronounced minima, the frequency of oscillations becomes lower and lower until the phase state is completely trapped 
in one of the two minima. 
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Figure Captions 

Figure 1: Energy spectrum for the Lipkin model with N =10 particles as a function of the strength parameter \- 
The positive part of the spectrum is obtained by reflecting the negative part around zero. 

Figure 2: The gap function A = E\ — Eq as a function of the strength parameter \- 

Figure 3: First three numerical derivatives of the gap curve presented in the previous figure. The vertical bars 
indicate the relevant values of \ where the correlations change the gap in a peculiar way. In this figure the solid line 
indicates the first derivative (n = 1), the dashed line corresponds to n = 2, and the long dashed-short dashed line 
curve to n = 3. 

Figure 4: Family of potential curves associated to different values of the strength parameter \- The curve for 
X = shows the potential for a spin system without two-body correlations. The curves for \ = 1-2 and \ = 1-8 show 
the potential for the particular values indicated in the previous figure. The curve for \ = 2-5 exhibits a typical case 
of spin tunneling. 

Figure 5: Family of curves showing the probability of finding the symmetric combination as a function of time for 
some values of the strength parameter x- We have prepared a symmetric combination at t = 0. As x increases the 
frequency diminishes. 

Figure 6: Family of curves showing the mean phase value as a function of time for some values of the strength 
parameter x for the symmetric combination. It is shown that, as x increases, the system state bounces with a lower 
frequency between the two wells. 
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